A statistical approach to the QCD phase transition 
— A mystery in the critical temperature 
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We study the QCD phase transition based on the statistical treatment with the bag-model picture 
of hadrons, and derive a phenomenological relation among the low-lying hadron masses, the hadron 
sizes and the critical temperature of the QCD phase transition. We apply this phenomenological 
relation to both full QCD and quenched QCD, and compare these results with the corresponding 
lattice QCD results. Whereas such a statistical approach works well in full QCD, it results in an 
extremely large estimate of the critical temperature in quenched QCD, which indicates a serious 
problem in understanding of the QCD phase transition. This large discrepancy traces back to 
the fact that enough number of glueballs are not yet thermally excited at the critical temperature 
T c ~ 280 MeV in quenched QCD due to the extremely small statistical factor as e" mG/Tc ~ 0.00207. 
This fact itself has a quite general nature independent of the particular choice of the effective model 
framework. We are thus arrive at a mystery, namely, what is really the trigger of the deconfinement 
phase transition. 
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I. INTRODUCTION 

The quark-gluon-plasma (QGP) is one of the most in- 
teresting targets in the finite-temperature quark-hadron 
physics Q. Currently, the QGP creation experiment is 
being performed in RHIC project at BNL, and much 
progress in understanding the finite-temperature QCD is 
desired. Historically, the instability of the hadron phase 
was first argued by Hagedorn before the discovery of 
QCD. He pointed out the possibility of a phase transi- 
tion at finite temperature, based on the string or the 
flux-tube picture of hadrons^, H|. After QCD was estab- 
lished as the fundamental theory of the strong interac- 
tion, this transition was recognized as the deconfinement 
phase transition to the QGP phase, where quarks and 
gluons are liberated with the restored chiral symmetry. 
The QCD phase transition has been studied using vari- 
ous infrared effective models of QCD such as the linear 
a model the Nambu-Jona-Lasinio model the dual 
Ginzburg-Landau theory^] and so on. 

In order to study nonperturbative features of the QCD 
phase transition, the lattice QCD Monte Carlo calcula- 
tion serves as a powerful tool directly based on QCD. It 
has been already extensively used to study the nature of 
the QCD phase transition. At the quenched level, SU(3) 
lattice QCD indicates the existence of the deconfinement 
phase transition of a weak first order at T c ~ 260 — 280 
MeV 0. On the other hand, in the presence of dynam- 
ical quarks, it indicates the chiral phase transition at 
T c = 173(3) MeV for N f = 2 and T c = 154(8) MeV for 
Nf = 3 in the chiral limit |8| . 

In this paper, we attempt to understand the physi- 



cal implications of the recent lattice QCD results on the 
QCD phase transition, and point out an abnormal nature 
of the quenched QCD phase transition based on the sta- 
tistical argument. In the actual calculation, we adopt the 
statistical approach to the QCD phase transition with 
the bag-model picture of hadrons M. Although it is a 
simply-minded phenomenological effective model, it can 
reproduce the critical temperature of the full-QCD phase 
transition "amazingly" well. In spite of this success, this 
approach terribly overestimates the critical temperature 
in quenched QCD. We consider the essential cause of this 
overestimate, and point out a serious problem hidden in 
the QCD phase transition in a model-independent man- 



The contents are organized as follows. In Sect. |J, 
we give a brief review of the statistical approach to the 
QCD phase transition with the bag-model picture of 
hadrons, and derive a phenomenological relation among 
the hadron masses, the hadron sizes and the critical tem- 



perature. We apply this relation to full QCD. In Sect. Ill 



we apply this relation to quenched QCD without dynam- 
ical quarks, and point out a serious problem on the crit- 
ical temperature. In Sect. IV, after the summary of the 



results, we attempt to clarify the essence of this prob- 
lem, and discuss an abnormal nature of the deconfine- 
ment phase transition in quenched QCD. 
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II. A STATISTICAL APPROACH TO THE 
FULL-QCD PHASE TRANSITION 

We investigate analytically the features of the QCD 
phase transition based on the statistical treatment with 
the bag-model picture of hadrons. We begin by deriving 
a phenomcnological relation between the critical temper- 
ature T c and the properties of the hadrons, i.e., the mass 
and the size. In the bag-model picture, quarks and gluons 
are assumed to be confined inside a spherical bag. Here, 
color confinement is simply taken into account through 
the bag- like intrinsic structure of hadrons. At low tem- 
perature, only a small number of such bags are ther- 
mally excited, and the thermodynamic properties of the 
system are described in terms of these spatially isolated 
bags With the increasing temperature, the number of the 
thermally excited bags increases. Gradually, these bags 
begin to overlap one another, and they finally cover the 
whole space region at a critical temperature T c in this 
picture. Above T c , since the volume of the outside dis- 
appears, the whole space is filled with liberated quarks 
and gluons, and the thermodynamic properties of the sys- 
tem is now governed by quarks and gluons. In this way, 
the QCD phase transition is described in terms of the 
overlaps of the thermally excited bags. Such a physical 
interpretation of the QCD phase transition is called as 
the closed packing picture. 

To proceed, we define the spatial occupation ratio 
rv(T) at temperature T to be the ratio of the total vol- 
ume of the spatial regions inside the thermally excited 
bags to the volume V of the whole space region. In 
the closed packing picture of the QCD phase transition, 
ry(T) plays the key role, which is estimated as 



(1) 



VA 7?^ f d 3 k 
^ n n 3 J (2nf 



) 3 e \A»'+£ 2 /T 



= ^A tl ^T 3 /K/T), 
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and its functional form is plotted against m = m/T in 
Fig. H Note that /(m/T) is a characteristic function to 
describe the thermal contribution of the boson with the 
mass m at the temperature T |ic[| . For m 3> T, /(m/T) 
decreases exponentially with m/T, and the thermal con- 
tribution is expected to become negligible. 

As mentioned before, in the closed packing picture, the 
phase transition takes place, when the thermally excited 




FIG. 1: The function f(m) in Eq. 
m/T. 
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plotted against m 



bags almost cover the whole space region. Hence, the 
critical temperature T c is estimated by solving 



rv(T c ) = 1. 



(3) 



Now, we investigate the phase transition in full QCD 
using the statistical approach. In full QCD, the lightest 
physical excitation is the pion, and all the other hadrons 
are rather heavy as m » m v , T c . In fact, the pion is 
considered to play the key role in describing the thermo- 
dynamic properties of full QCD below T c from the view- 
point of the statistical physics. Hence, in most cases, only 
the pionic degrees of freedom are taken into account in 
the hadron phase in the argument of the full-QCD phase 
transition. By using the isospin degeneracy A,r = 3, the 
mass TOtt = 140 MeV and the radius ~ 1 fm, we solve 
Eq. (§) with r v {T) = 3RlT 3 f(m^/T) to estimate the 
critical temperature as T c ~ 183 MeV. Considering its 
closeness to the full lattice QCD result with Nf = 2, i.e., 
T c ~ 170 MeV, the statistical approach to the full-QCD 
phase transition seems to be rather good. 

We now consider the m^-dependence of the critical 
temperature T c . Note that, in the actual lattice QCD 
calculations, the pion mass is taken to be still rather 
heavy as > 400 MeV for the technical reasons. From 
these data, the critical temperature T c in the chiral limit 
is obtained using the chiral extrapolation. In Ref.[g, the 
authors parametrized the full lattice QCD result of T c as 



(TTt \ 



where a denotes the string tension. Strictly speaking, the 
phase transition becomes just a cross-over for intermedi- 
ate values of m^. Hence, in Ref.||, the pseudo critical 
temperature is adopted as T c , which is determined from 
the peak positions of the susceptibilities of the Polyakov 
loop and so on. 
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FIG. 2: The critical temperature T c plotted against the pion 
mass m-n in the Nf — 2 case. The dashed curve denotes the 
result retaining only the contribution from pions. The solid 
curve denotes the result including also contributions from p 
and uj mesons. The triangle denotes the lattice data taken 
from Ref. jTl| obtained with standard (std) fermion action. 
The square denotes the lattice data taken from Ref. || ob- 
tained with the improved staggered (p4) fermion action. 



FIG. 3: The critical temperature T c plotted against the 
pseudo-scalar meson mass tops in the SU(3) / symmetric case. 
The dashed curve denotes the result retaining only the con- 
tribution of the pseudo-scalar mesons. The dot-dashed and 
the solid curves denote the results including also the contri- 
bution of the flavor-octet vector mesons with R = 1 fm and 
R = 0.75 fm, respectively. The circle denotes the lattice data 
taken from Ref.[g] obtained with the improved staggered (p4) 
fermion action. 



For the Nf = 2 case, we plot, in Fig. the theoretical 
estimate of T c against based on Eq. (^) in the statis- 
tical approach with the bag-model picture. The dashed 
curve denotes T c retaining only the contribution from 
pions with R^ ~ 1 fm. The solid curve denotes the re- 
sult in the case including also the contributions from the 
low-lying vector mesons, such as p and u>, which are the 
next lightest particles in Nf = 2 full QCD. Here, we 
have used A p = 3 x 3 = 9, m p = 770 MeV, R p ~ 1 
fm, A w — 3, m u = 783 MeV, R u ~ 1 fm as inputs, 
which are treated as m^-independent constants. These 
vector mesons give an additional contribution to ry{T) 
as Sr v (T) = 9R 3 p T 3 f(m p /T) + 3^T 3 f(m u /T). The 
triangle and the square in Fig. ^ denote the lattice data 
taken from Refs.[||, pj| . In spite of its roughness, the re- 
sulting behavior of the critical temperature T c seems to 
coincide rather well with the lattice QCD result. In this 
calculation, the low-lying vector mesons, p and uj, give a 
small but non-negligible contribution to the critical tem- 
perature. 

We next investigate the idealized SU(3)/ symmetric 
case to analyze the other lattice QCD data. In this 
case, the pseudo-scalar (PS) octet mesons, such as pions, 
kaons and ryg, are the lightest and possess the same mass 
tops in common with the degeneracy as Apg = 8. We 
plot, in Fig. the critical temperature T c against mpg. 
The dashed curve denotes T c retaining only the contri- 
bution from PS-mesons with i?p§ ~ 1 fm, which leads 
to r v (T) = 8Rl s T 3 f(m PS /T). The dot-dashed curve 
in Fig. H denotes the results including also the contribu- 
tions from the octet vector mesons with Ay = 8 x 3 = 24, 
my ~ 770 MeV, Ry ~ 1 fm. (Inclusion of the flavor- 
singlet vector meson does not change the result so much.) 



These vector mesons give an additional contribution to 
rv(T) as Sr v (T) = 2AR Y T Z f(m v /T). The circle in 
Fig. H denotes the lattice data taken from Ref.(|]. We 
see that both the dashed and the dot-dashed curves are 
roughly consistent with the lattice QCD results. We note 
that the small deviation almost disappears by slightly ad- 
justing the bag size as i?ps = Rv = 0.75 fm, as is shown 
in Fig. H with the solid curve. 

In this way, this simple statistical approach with the 
bag-model picture works rather well in reproducing the 
critical temperature T c of the full-QCD phase transition. 
Here, the main contribution is given by the thermal pion, 
and the low-lying vector mesons give a small but non- 
negligible contribution to T c . 



III. A MYSTERY IN THE QUENCHED QCD 
PHASE TRANSITION 

In spite of the absence of the dynamical quarks, 
quenched QCD provides various important nonpertur- 
bative features such as the color confinement and instan- 
tons. Hence, to understand some of the nonperturbative 
natures of QCD, quenched QCD plays the primary role, 
serving as a simplified version of the complicated real 
problems. In this section, we apply the statistical ap- 
proach with the bag-model picture to quenched QCD in 
a similar manner done in the previous section. 

Due to the color confinement, only the color-singlet 
modes can appear as physical excitations, and such 
modes are called as glueballs in quenched QCD. The mass 
spectrum of the glueballs is known through the quenched 
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FIG. 4: The critical temperature T c plotted against the 
lowest ++ glueball mass mofs) in the statistical approach 
assuming various glueball sizes as Rq(S) = 0.5, 1.0, 1.5 fm. 
The cross (x) indicates the quenched lattice QCD results, 
T c = 280 MeV and m G(s) = 1730 MeV. 



lattice QCD calculations [||, [b| 0. The lightest phys- 
ical excitation is the ++ glueball with m G = 1.5 — 1.7 
GeV. Following the similar argument as in the full QCD 
case, the lowest ++ glueball is expected to play the 
key role in describing the thermodynamic properties of 
quenched QCD in the confinement phase. Hence, we first 
take into account only the ++ glueball. As for the size 
R of the scalar glueball, we note that there is no widely 
agreed value on it. We adopt a rather small value as R ~ 
0.4fm, which has been indicated by recent lattice QCD 
studies 0,0. 

We use the degeneracy Aq(s) = 1) the mass ttig(s) — 
1730 MeV @, and R G{S) ~ 0.4 fm as inputs. Then, 
the spatial occupation ratio is given by ry{T) = 
^G(S) ^ f( m G(S)/'^')j an( i the critical temperature is es- 
timated as T c ~ 827 MeV from Eq. (§). This estimate 
is too much larger than the quenched lattice QCD re- 
sult as T c = 0.629(3)^ ^ 280 MeV in Ref.g with 
= 450 MeV. (In other words, only a tiny fraction 
of the space region is covered by the thermally excited 
bags as r v (T) = 0.0021 at T = 280 MeV.) This large 
discrepancy would be a serious problem in the quenched 
QCD phase transition. 

To seek a possible solution, we examine the case with 
a larger value of the glueball size R, since there is no 
established value on R. In Fig. [|, the critical temperature 
T c is plotted against the lightest ++ glueball mass toq(s) 
in the statistical approach for various glueball sizes as 
^G(S) = 0.5,1.0,1.5 fm. The cross (x) indicates the 
quenched lattice QCD results, T c — 280 MeV and to g = 
1730 MeV. In this argument, to reproduce the critical 
temperature as T c ~ 280 MeV, the glueball size is to be 
abnormally large as Rg(s) — 3-1 fm in the vicinity of T c . 
However, such a drastic thermal swelling of the lowest 
scalar glueball was rejected by the recent lattice QCD 



studies |l5|, [l(|, which states that the thermal glueball 
size is almost unchanged even near T c . In any case, the 
glueball size does not seem to provide the solution on this 
discrepancy. 

We may seek for the solution in the drastic pole-mass 
reduction of the ++ glueball near the critical temper- 
ature as was suggested in Ref.(6) in the context of the 
dual Ginzburg-Landau theory. In Fig. ^, we see that, for 
the problem to be settled, the pole-mass reduction must 
be as significant as m G (s)(Tc) < 500 MeV. However, in 
the recent lattice QCD calculations |ll| |6), it has been 
reported that the thermal ++ glueball persists to hold 
a rather large pole-mass as m G (g)(T ~ T c ) ~ 1250 MeV. 
Hence, we have to seek for another possibility by includ- 
ing the contribution of the excited-state glueballs. 

In addition to the lowest ++ glueball with me(S) — 
1730 MeV and Rg(s) — 0.4 fm, we consider the ther- 
mal contribution from the lowest 2 ++ glueball, which 
is the next lightest hadron in quenched QCD. We take 
A G(T) = 5, m G(T) ~ 2400 MeV [jfj, R G{T) ~ 1 fm. 
(Here, we assume it to have a typical hadron size.) 
Then, the spatial occupation ratio receives a correction 
as Srv(T) = 5i? G ^ T jT 3 /(m G ( T )/T), and the correction 
amounts to Sr v (T) = 0.0223 at T = 280 MeV. The re- 
sulting critical temperature is given as T c ~ 432 MeV, 
which is still too large. We note that the realistic glue- 
ball size would be more compact. However, if so, its 
contribution becomes more negligible. 

Besides these two low- lying glueballs, the following ex- 
cited states are predicted in Ref.jjjj as h (2590), 
0*++(2670), l+-(2940), 2-+(3100), 3+^(3550), 
0*-+(3640), 3++(3690), 1~ (3850), 2*-+(3890), 
2— (3930), 3— (4130), 2+-(4140), 0+~(4740). We 
include all the contribution from these excited states, 
assuming the unknown glueball size as a typical hadron 
size, i.e., Rq ~ 1 fm. The correction by these excited 
states amounts to only Sr v (T) = 0.0113 at T = 280 
MeV, and the resulting critical temperature is estimated 
as T c = 395 MeV, which is again too large. 

In this way, we observe that the statistical approach 
does not work at all in quenched QCD. The direct cause 
of this failure is the extremely small statistical factor 
e -™G(s)/ T c ro 0.0021, which strongly suppresses the exci- 
tations of the glueballs. As a consequence, only a remark- 
ably tiny fraction of the space region can be covered by 
the thermally excited bags of glueballs at T = 280 MeV, 
which cannot be the driving force of the QCD phase tran- 
sition. 



IV. SUMMARY AND DISCUSSIONS —WHAT 
IS THE TRIGGER OR THE DRIVING FORCE 
OF THE QCD PHASE TRANSITION ? 

We have analytically studied the QCD phase transition 
based on the statistical treatment with the bag-model 
picture. We have derived a phenomenological relation 
among the critical temperature T c , the mass and the size 
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of the low- lying hadrons. First, we have applied this 
relation to full QCD, and have compared the analytical 
results with the full lattice QCD results. In full QCD, 
we have found that this approach works amazingly well. 
As is expected, the pionic contribution seems to play the 
main role in determining the critical temperature, and 
the low-lying vector mesons are found to provide a small 
but non-negligible contribution. 

Unlike full QCD, we have found that the statistical ap- 
proach terribly overestimates the critical temperature as 
T c ~ 827 MeV in quenched QCD, and that only a tiny 
fraction of the space region is covered by the thermally 
excited bags of glueballs at T ~ 280MeV. We have con- 
sidered the possibility of a swelling of the glueball size 
and reduction of the glueball mass near the critical tem- 
perature. However, both of these two have not provided 
us with the solution on the large discrepancy. In spite of 
including all the contributions from the 15 low-lying glue- 
balls up to 5 GeV predicted in quenched lattice QCD jf|] , 
the discrepancy remains to be still large. In other words, 
the number of the thermally-excited glueballs is still too 
small at T ~ 280MeV, even after so many glueball ex- 
cited states are taken into account. The direct origin of 
this discrepancy is the strong suppression of the thermal 
excitation of glueballs due to the extremely small statisti- 
cal factor as e ~ m G(s)/ T = 0.00207 at T = 280 MeV even 
for the lightest glueball, which leads to the insufficient 
amount of the covered space by the thermally excited 
bags of glueballs. 

Although several arguments given so far have been 
based on the bag-model picture, this problem itself has 
a quite general nature, which can go beyond the reliabil- 
ity of the model framework. We finally reformulate this 
problem in a model-independent general manner by con- 



sidering the inter-particle distance I instead of the bag 
size R. Note that the reliability of the statistical argu- 
ment becomes improved in the dilute glueball gas limit. 
Now, taking into account all the low-lying 15 glueball 
modes up to 5 GeV predicted in quenched lattice QCD 
p2[ , we calculate the inter-particle distance I of the glue- 
balls based on only the statistical argument. At T = 280 
MeV, the inter-particle distance is estimated as I ~ 5 
fm. It follows that the deconfinement phase transition 
takes place at T ~ 280MeV, when the glueball density 
becomes p = l/{4f (2.5fm) 3 } ~ l/(4.0fm) 3 . Remem- 
ber that the theoretical estimates of the glueball size are 
rather small as i? G (s) < 0.4 fm |l5|, p|, |P7|]. Hence, this 
density seems to be rather dilute. Since the long-range 
interaction among glueballs is mediated by the virtual 
one-glueball exchange process in quenched QCD, the in- 
teractions among glueballs are exponentially suppressed 
beyond its Compton length l/mG(S) — 0.112 fm. There- 
fore, one cannot expect the strong long-range interaction 
acting among the spatially-separated thermal glueballs. 

In this way, it is quite difficult to imagine how such 
a too rare excitation of thermal glueballs can lead to 
the phase transition. We are thus arrive at the mys- 
tery What is really the trigger or the driving force of 
the deconfinement phase transition in quenched QCD ? 
In order to understand the QCD phase transition, this 
problem should be seriously considered. 
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